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ABSTRACT 

We use Hubble Space Telescope/Space Telescope Imaging Spectrograph 
(HST/STIS) longslit low-resolution spectroscopy from 1150 A to 10,300 A to 
study the physical conditions in the narrow-line region (NLR) of the Seyfert 2 
galaxy Markarian 3. We find from the He II A1640/A4686 line ratio and the 
Balmer decrement that the extinction within Markarian 3 along the line-of-sight 
to the NLR is best characterized by a Large Magellanic Cloud (LMC) type ex- 
tinction curve. We observe an extinction gradient increasing from west to east 
along the STIS slit (at position angle 71° measured east from north) in both 
line and continuum emission. We infer from this gradient that the host galaxy 
disk is tilted towards the observer in the east: the line-of-sight to the eastern 
emission-line cone intersects more dust in the plane of the galaxy than that to 
the western cone. From emission-line diagnostics we find that the NLR gas is 
photoionized by the hidden active galactic nucleus (AGN) continuum and that 
its density decreases with increasing distance from the center. We model the 
observed continuum as a combination of reddened host galaxy light from an old 
stellar population, reddened H + and He +2 recombination continua, and less red- 
dened scattered light from the central engine with spectral index a—1 (L„ oc v~ a ). 
The host galaxy to scattered-light ratio is estimated to be 3:1 at 8125 A in O'.'l 
x l'/8 aperture. Using a two-component power-law model for the ionizing con- 
tinuum (a=2 for 13.6 eV < E < 0.2 keV and a=l for 0.2 keV < E < 50 keV) 
we find that the covering factor (normalized for our observation aperture) of the 
NLR gas is >0.7%. We estimate that the amount of intrinsic non-ionizing UV 
continuum scattered into our line-of-sight is 0.04%. This is consistent with our 
estimate of the scattering fraction for broad C IV AA1548,1551 emission. 
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1. Introduction 

Markarian 3 is classified as a Hubble-type SBO galaxy (Adams f977) with a Seyfert 2 
(Sy 2) active galactic nucleus (AGN) (Weedman & Khachickian 1969). Its spectrum is 
characterized by narrow forbidden and permitted emission lines and excess UV flux compared 
to normal galaxies (Weedman 1973). Thompson et al. (1980), Schmidt & Miller (1985), 
Miller & Goodrich (1990) and Tran (1995a) observed broad permitted lines and AGN 
continuum emission in its polarized spectrum. These observations are part of the foundation 
for the "unified" model for active galaxies (Antonucci 1993; Urry & Padovani 1995). In 
this model the presence or absence of broad permitted-emission lines and a bright, variable 
AGN continuum in the observed (unpolarized) spectrum depends on the viewing angle to the 
central engine. The systemic velocity is 4050 km s" 1 (z=0.0135) based on H I 21 cm emission 
(Tifft & Cocke 1988). The distance to Markarian is 53 Mpc if H =75 km s" 1 Mpc -1 . At 
this distance, 1" corresponds to 257 pc. 

Markarian 3 has been well studied from X-rays to radio wavelengths. Its NLR has a 
reverse- "S" shape greater than 2" in extent (Capetti et al. 1995). In a study of near-UV to 
far-IR absorption lines, Gonzalez Delgado, Heckman & Leitherer (2001) found evidence of 
an old stellar population in the nucleus of Markarian 3. A 2"-long radio jet detected in 5- 
GHz MERLIN observations (Kukula et al. 1993) has a similar shape and orientation to the 
NLR morphology determined by HST/Faint Object Camera (FOC) narrow-band imaging 
observations centered on the [O II] A 3727, [O III] A 5007, Ha, and H7 emission-lines (Capetti 
et al. 1995). Its hard X-ray spectrum can be modeled as Compton reflection from molec- 
ular gas in a torus that obscures the central source of ionizing radiation (Georgantopoulos 
1999; Turner et al 1997). However, Georgantopoulos et al. (1999) found a factor of two 
variation in the X-ray source flux between 4 keV and 20 keV over a 200-day period in RXTE 
observations. They inferred that a fraction of the intrinsic X-ray continuum penetrates the 
putative obscuring torus. A comparison of ASCA, GINGA, and BBXRT data revealed that 
the extended soft X-ray flux (E<4 keV), which covers the same ~12" range as the extended 
narrow-line- region (ENLR) (Sako et al. 2000), has shown no such variation over ~13 years 
(Iwasawa et al. 1994). 

Ruiz et al. (2001) studied the kinematics of the Markarian 3 NLR using [O III] A5007 
measurements from a STIS/CCD slitless spectrum and from the same longslit data that is 
presented in this paper. We use the results of their kinematic analysis in this paper. The 
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extracted spectra from several locations along the longslit show evidence of two kinematic 
components: one is redshifted and the other is blueshifted in the rest-frame of Markarian 3. 
These components are interpreted as line-of-sight observations of radially outflowing gas on 
opposite sides of a NLR bi-cone. Their best-fit kinematic model for a NLR with an angular 
extent of 2" and inner and outer opening half-angles of 15° and 25°, respectively, is radial 
acceleration from the nucleus to a turnover distance (measured along the bi-cone axis) at 
~80 pc followed by deceleration to the systemic velocity at ~260 pc. The maximum space 
velocity in the model is 1750 km/s. In this model the bi-cone is tilted towards the observer 
in the east and away from the observer in the west by 5 degrees. The position angle of the 
the bi-cone is 70° east of north (Schmitt & Kinney 2000). We note that the host galaxy 
disk major-axis position angle is 28° and the inclination of the disk is 33° (Schmitt & Kinney 
2000). We determine the direction of the host-galaxy tilt in this paper. 

In this paper we present an analysis of longslit spectra obtained in each of the HST/STIS 
low-resolution modes: G140L, G230L, G430L, and G750L. These spectra show fifty emission- 
lines from ~21 ions. The brightest emission-lines have an angular extent of ~3'/4 or a 
projected linear extent of 870 pc. In §2 we describe the observations and data reduction. In §3 
we describe the emission-line measurement procedure and discuss the extinction corrections. 
In §4 we present temperature, density and ionization diagnostics. In §5 we show the observed 
continuum and describe its constituents. We also estimate the covering factor of the NLR 
gas, and the scattering fraction for the continuum produced by the hidden AGN. Our physical 
picture of the NLR, including the relative orientation of the host galaxy disk, is presented 
in §6. We summarize our results in §7. In the forthcoming paper II, we will use these 
observations and the photoionization model code CLOUDY by Ferland (1988) to constrain 
the sizes, densities, and spatial distribution of gas clouds in the NLR. 



2. Observations and Data Reduction 

HST/STIS longslit spectra were obtained on 22 August 2000 in all four low- resolution 
modes, G140L, G230L, G430L, and G750L providing wavelength coverage from 1150 A to 
10,300 A. The observations are listed in Table 1. The 52"xO / .'l slit was used for all spectral 
modes. This aperture is superposed on the target acquisition image in Figure 1. The slit 
was oriented to position-angle 71° east of north to match the elongated reverse- "S" shaped 
feature in the NLR. East is towards the bottom of this figure, and North is towards the 
left, as indicated by the arrows. The target acquisition image was obtained through the 
STIS long-pass F28X50LP filter, which spans red to near-IR wavelengths. The bright knots 
in the image are primarily emission from lines that fall within the wavelength range of the 
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F28X50LP filter. The brightest lines in the visual and near-IR are [O I] AA6300, 6363, 
[N II] AA6548, 6584, Ha A6563, [S II] AA6716, 6731, and [S III] AA9069, 9532. 

The spectra were processed with the CALSTIS IDL software developed by Lindler 
(1998) for the STIS Instrument Definition Team. This software is used to perform bias and 
dark subtraction, hot pixel correction, flat-field division, and for the CCD images, cosmic-ray 
rejection. Hot pixels were identified in STIS dark images and corrected by interpolating in 
the dispersion direction. Cosmic-ray artifacts were removed in the CCD data by comparing 
the multiple image readouts obtained within each spectral mode. The multiple readouts 
obtained in the G140L, G430L, and G750L modes were combined to increase their signal- 
to-noise ratios. The wavelength scales were corrected for zero-point shifts using hole-in- 
the-mirror (HITM) lamp calibration exposures made after each science exposure. The two 
G140L spectra have a relative wavelength offset of 0.15 A , or 0.25 pixel. Both spectra were 
shifted to match a central wavelength of 1418.821 A. No such shifts were required for the 
mode G430L and mode G750L input spectra. 

The spectra were processed in "extended-source" mode. In this mode, each column of 
the output spectral image corresponds to a unique wavelength and each row corresponds to 
a unique spatial position. This rectification was performed using a bi-linear interpolation al- 
gorithm (Lindler 1998). The units of the output spectral image are ergs s _1 cm -2 A -1 cross- 
dispersion-pixel -1 . In addition to this rectification, all spectral images were transformed to 
the same spatial plate scale (1 pixel = 0'.'05071) and alignment as that of spectral mode 
G430L. This simplifies comparison of the same NLR clouds observed in different spectral 
modes. The co-alignment among modes in the spatial direction was determined by cross- 
correlating the slit fiducial marks that are evident in the respective wavelength calibration 
spectra. 

Finally, a background correction was made for each of the spectral images. For modes 
G230L, G430L, and G750L, the background subtraction algorithm follows that described 
by Lindler (1998) for point-source extractions. Briefly, two eleven-pixel wide strips are 
selected a suitable distance away (typically more than 200 pixels, or 10") from the center 
of the NLR spectrum. Each strip is averaged in the spatial direction, then the two average 
background spectra are in turn averaged together. The resultant background spectrum was 
smoothed three times: once with a median filter, and again with two different mean filters. 
This smoothed background spectrum was subtracted from each line of the spectral image. 

The background subtraction for mode G140L required more care due to contamination 
from geo-coronal Lyct A1216 and O I A1303 emission . The pixels in the averaged, smoothed 
background spectrum that contain the two geo-coronal lines were replaced with pixels from 
the averaged, unsmoothed background spectrum. This method minimizes subtraction resid- 
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uals at the positions of these lines. 

We examined the impact of fringing on the CCD/G750L mode spectrum. Fringing in 
the STIS/CCD G750L mode spectra can reach peak-to-peak amplitude variations of greater 
than 20% at A > 9000 A (Proffitt et al. 2002; Malumuth et al. 2003). The fringing is 
caused by interference among reflected waves in the CCD substrate. We investigated the 
fringe amplitude of our data using a tungsten lamp flat-field (Plait 2000) and a model 
flat-field (Malumuth et al. 2003). We concluded that the uncertainty in the data in the 
wavelength region affected by fringing (A > 7500 A) is dominated by Poisson noise rather 
than fringing. No correction for fringing was deemed necessary. 

We display l'/8 sub-image extractions from the co-aligned, background-subtracted spec- 
tral images for G140L, G230L, G430L and G750L in Figures 2, 3, 4, and 5, respectively. 
The overlaid 1-D spectrum is the sum of the spectral image in the cross-dispersion direction. 
The spatial orientation of each spectral image is West at the top of the figure and East at 
the bottom of the figure, consistent with Figure 1. We detected over 60 emission lines from 
ionic species that require ionization potentials up to 100 eV. 

3. Emission-Line Measurements 

3.1. Integrated Fluxes 

We made seven discrete measurements of each emission line in the spatial direction. 
Individual spectra were extracted by summing array rows in selected bins in the spatial 
direction. The bin sizes were selected to extract spectral information for individual clouds 
and to achieve a reasonable signal-to- noise ratio for the the He II A1640/A4686 measurement 
which is used to determine the extinction along the line-of-sight to the NLR (see §3.2). 
The bin sizes range from 0.2-0.3 arcseconds, or 52-77 pc, using the adopted distance scale 
described in the Introduction. The full spatial range along the slit spanned by the defined bins 
was also constrained by the spatial extent of the two He II emission lines. Several emission- 
lines have spatial extents greater than the analyzed l'/8 range, but the extinction correction 
outside this range cannot be determined from the low signal-to-noise He II emission. 

The measurement bins are shown in Figure 6. The bins were selected by visually com- 
paring the spatial profile of each emission-line with that of the target acquisition direct 
image. The arrow on each row of displayed emission lines indicates bin O'.'O. It contains the 
kinematic zero-point in the Ruiz et al. (2001) velocity map. This coincides with the posi- 
tion of the hidden-nucleus derived by Kishimoto et al. (2002) from UV-imaging polarimetry 
data. 
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The line flux measurement procedure began by removing all lines except the line of inter- 
est from the one-dimensional spectral extraction. We executed this removal by interpolating 
between continuum values adjacent to each unwanted line. We then performed a linear fit to 
the remaining continuum. We subtracted this fit from the spectrum, then fitted a Gaussian 
function to the line of interest. The uncertainties in the line fluxes were estimated by making 
three such measurements for each line, each time interactively selecting different continuum 
regions. The wavelength range of the continuum regions varied between the full spectral 
range of the spectral mode of interest (G140L: 611 A, G230L: 1616 A, G430L: 2800 A, and 
G750L: 5030 A) and one-third of the the full range. The reported line flux is the average 
value of the three measurements. The line flux uncertainty is the standard deviation of the 
three measurements added in quadrature with the uncertainties from Poisson statistics. This 
uncertainty is dominated by the standard deviation of the three line flux measurements. 

In any measurement bin where two kinematic components and/or two different lines 
are blended together, we used the [O III] A5007 solution for that bin as a template. There 
were no compelling incidences where three or more components were required to fit the mea- 
surements using the selected bins. We assumed that the kinematics of the [O III] A5007 gas 
matches that of the line of interest. The [O III] A5007 measurements with two components 
were fitted with six free parameters, or three for each of the two components: the peak 
flux, the Gaussian standard deviation, and the central wavelength. The Gaussian standard 
deviations were converted to velocity dispersions. The fitted central wavelengths were used 
with the rest wavelength to compute the radial velocities of the components. The dispersions 
were converted to Gaussian standard deviations for the line(s) of interest. These standard 
deviations were corrected to match the observed line profiles by adding in quadrature the 
appropriate spectral resolution term from Table 1. The radial velocities were used to de- 
rive the wavelength offsets for the line(s) of interest. This procedure leaves only the peak 
fluxes of the components as free parameters in the Gaussian fits for the line(s) of interest. 
The lines deblended in this manner were [Ne III] A3869 from H8 A3889, [S II] A4074 from 
R5 A4102, H 7 A4340 from [O III] A4363, [Fe VII] A5159 from [N I] A5200, [Fe VII] A5722 
from [N II] A5755, [O I] A6300 from [O I] A6363, and He I A7065 from [Ar III] A7136. For the 
Ha A6563+[N II] AA6548,6584 blend we fixed the [N II] A6584/A6548 flux ratio to 3:1 in ad- 
dition to the aforementioned constraints. The 3:1 ratio applies to a ~10,000 K nebulae with 
density less than 10 5 cm -3 . This is the critical density required to collisionally de-excite the 
upper level ( 1 D2) responsible for the [N II] AA6548,6584 doublet emission. We used the same 
procedure to deblend the C IV AA1548,1551 doublet and the Mg II AA2796,2803 doublet. 
For both of these doublets the short-to-long wavelength line-flux ratio constraints are 2:1 for 
optically thin lines. The spectral resolution is too low to separate the individual components 
of the remaining blended lines: Si IV A1398 multiplet with the O IV] A1402 multiplet, and 
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[S II] A6716 with [S II] A6731. Figure 7 shows the fit for the two- component [O III] A5007 
line in the central bin. We list flux measurements relative to H/3 A4861 for forty-six emission 
lines at seven spatial positions in Table 2. The la uncertainties in the observed flux ratios 
are listed in parentheses below each measurement. 

3.2. Extinction Correction 

We corrected the observed flux measurements for both Galactic extinction along the 
line-of-sight to Markarian 3, and intrinsic extinction caused by dust within Markarian 3. 
The Galactic color-excess of E(B-V) = 0.19 was obtained from the dust map of Schlegel et 
al. (1998). We used the Savage & Mathis (1979) extinction curve to correct for Galactic 
extinction. 

We used the He II A1640/A4686 line ratio (Seaton 1978) to determine the extinction 
within Markarian 3 along the line-of-sight. This ratio has a longer wavelength baseline than 
the Balmer decrement (Ha A6563/H/3 A4861) and therefore gives a better estimate of the 
reddening. It is also insensitive to collisional excitation. The ratio can vary from 5.6 for a 
low density (n e — > cm" 3 ) 5,000 K plasma to 7.6 for a high density (n e ~10 6 cm" 3 ) 20,000 K 
plasma (Seaton 1978). We used a ratio of 7.2. This is appropriate for a 10,000 K nebula 
with Case B recombination and n e < 10 4 cm -3 (Seaton 1978). We calculated the color 
excess for each measured NLR component using the following screen geometry relation 

E(B -V) = 2.5 x log w (7.2 x (/ 4686 //i64o))/(*i640 - X 4686 ) (1) 

where f 4686 and fi 640 are the Galactic extinction corrected He II line fluxes and X 1640 and 
X 4686 are the normalized extinction values (X A = E(A - V) / E(B-V)) at the wavelengths 
of the He II lines. The la uncertainties in the color excess values are on the order of the 
range of color excess values calculated using the theoretical upper and lower limits of the 
He II A1640/A4686 line ratio described above. 

To determine the best extinction curve for the NLR of Markarian 3, we calculated E(B- 
V) using several different curves from the literature, applied the different color excess values 
to the Ha A6563 and H[3 A4861 flux measurements, and compared the resultant Balmer 
decrement values with theoretical values. The Balmer decrement can vary from 2.74 for a 
high density (n e = 10 6 cm" 3 ) 20,000 K gas to 3.05 for a low density (n e — > 0) 5,000 K gas 
assuming Case B recombination (Osterbrock 1989). We also required physically plausible 
values (typically <36 [Osterbrock 1989]) for the extinction corrected Lya A 12 16 to R(3 A4861 
flux ratios. We used following extinction curves: Galactic (Savage & Mathis 1979) , 6 2 Ori B 
(Bohlin & Savage 1981), the starburst curves of Calzetti et al. (1994) and Calzetti (1997), 
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SMC-bar (Witt & Gordon 1999), SMC (Hutchings 1982), average LMC (Nandy et al. 
1981), average LMC (Koornneef & Code 1981), and separate curves for within and outside 
of 30 Doradus in the LMC (Fitzpatrick 1985). 

The 6 2 Ori B, Galactic, starburst, SMC-bar curves yielded values too low to be phys- 
ically plausible. The Balmer decrement weighted mean values and their uncertainties are, 
respectively, 0.44±0.08, 2.33±0.09, 2.30±0.13, 2.36±0.10, and 2.02±0.11. The SMC curve 
(Hutchings 1982) yielded a Balmer decrement weighted mean value of 2.84±0.06, but the 
extinction corrected Lya A1216 values were higher than the theoretical upper limit by fac- 
tors of up to 2.8. The Nandy et al. (1981) and Fitzpatrick (1985) LMC curves (within 
and without 30 Doradus) produced weighted mean Balmer decrement values closer to the 
theoretical limits: 2.53±0.08, 2.52±0.08 (30 Dor), and 2.35±0.08 (non-30 Dor), respectively. 
We found that the LMC curve of Koornneef & Code (1981) gave a weighted mean value 
of the Balmer decrement (2.59±0.07) closest to the theoretical values and yielded physically 
plausible values for the Lja A1216 to H./3 A4861 flux ratios. We used this curve to de-redden 
the observed fluxes using a screen geometry in this paper. The extinction corrected line 
fluxes relative to H[3 A4861 are listed in Table 3. The la line-ratio uncertainties include a 
component due to the uncertainty in the extinction correction. 

The simple grain model of Mathis et al. (1977) for standard Galactic extinction 
curves attributes the continuous extinction to a size-distribution of various grains, and the 
2200 A bump to smaller (mainly graphite) grains. The empirically derived Koornneef & 
Code (1981) LMC curve exhibits a less pronounced 2200 A bump and increases more 
steeply with decreasing wavelength for A < 2200 A. The steeper rise of the Koornneef & 
Code (1981) LMC curve is attributed to an overall smaller distribution of dust grains than 
are required to produce the Galactic, starburst, and SMC-bar extinction curves described 
in this section. The agreement of the He II A1640/A4686 ratios and the Balmer decrement 
values with their theoretical values after applying the Koornneef & Code (1981) LMC curve 
suggests that the dust grain size distribution in the observed area of the Markarian 3 NLR is 
smaller than the standard Galactic distribution. Pitman & Clayton (2000) note that many 
active galaxies have dust extinction characterized by the absence of the 2200 A feature. The 
reddening curve is generally similar to that of the narrow-line Sy 1 Arakelian 564 (Crenshaw 
et al. 2002), but not as steep as that of the SMC (Hutchings 1982). 

The color excess is shown as a function of position in the NLR in Figure 8. There is 
a striking gradient in extinction increasing from west to east. This result might also have 
been inferred from inspection of the spectral images. Compare the spatial morphology of 
UV emission lines in Figures 2 and 3 with that of optical lines in Figures 4 and 5. Both 
redshifted and blueshifted points east of the nucleus are reddened, suggesting that the source 
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of extinction may be an intervening dust lane or screen between the observer and the NLR 
along the line-of-sight. The extinction gradient suggests that there is more dust in the screen 
along the line-of-sight to the eastern side of the NLR bi-cone than that of the western side, 
as observed along the STIS slit. Elmegreen & Block (1999) noted that dust in an inclined 
galactic plane can produce a similar extinction gradient in which the side tilted towards the 
observer is more reddened than the side that is tilted away (see further discussion on the host 
galaxy orientation in §6). We note that the magnitude of the extinction gradient is steeper 
when derived using an extinction curve that increases slowly with shorter wavelength in the 
UV, such as the 6 2 Ori B (Bohlin & Savage 1981) curve. The Hutchings (1982) SMC curve, 
which increases faster with decreasing wavelength in the UV, produces a shallower extinction 
gradient than that derived with the LMC Koornneef & Code (1981) curve. 

Kraemer et al. (2000) found no clear trend over their measured region (±3", or ~ 
±200 pcfor z=0.0033 and H =75 km s" 1 Mpc^ 1 ) of the NLR in NGC 4151 (Sy 1.5), although 
they note that the extinction may increase from the north-east to the southwest over the 
inner ±1" (~ ±60pc). They concluded that the lack of uniformity in extinction over the 
entire NLR suggests that the reddening is associated with the emission-line knots, rather 
than an external screen. Note that the orientation of the NGC 4151 NLR is such that all of 
the north-east side is redshifted and all of the south-west side is blueshifted (Nelson et al. 
2000). This is not surprising since the NLR bi-cone axis is tilted further toward our line-of- 
sight (35° out of the plane of the sky [Evans et al. 1993]) than that in Markarian 3. The 
host NGC 4151 galaxy is closer to face-on orientation than Markarian 3 with an inclination 
of 21° such that the eastern side is tilted towards the observer (Kinney et al. 2000). 

In NGC 1068 (Sy 2), Kraemer & Crenshaw (2000b) found in general that measurements 
of blueshifted clouds exhibited greater extinction than redshifted clouds. They inferred 
that both red and blueshifted components were observed through a foreground dust-screen, 
and that additional dust may be associated with the individual blueshifted clouds. The 
orientation of the NGC 1068 NLR is similar to that of Markarian 3, where blue- and redshifted 
clouds are observed along the same line-of-sight and the bi-cone axis is tilted towards the 
observer in the north-east by 5° (Crenshaw & Kraemer 2000). The angular extent measured 
along the the NLR bi-cone axis was ±5", or spatial extent ±360 pc (for z=0.0038 and 
H =75 km s _1 Mpc -1 ). The host galaxy disk is tilted towards the observer in the south 
by 28° (Kinney et al. 2000). The covering factors of the blueshifted clouds are expected 
to be small so that those clouds do not eclipse (and further redden) the redshifted clouds; 
therefore, the dust responsible for the observed extinction gradient in the Markarian 3 NLR 
is likely a screen outside the NLR. 

We note that most of the emission line profiles in the dual- velocity-component mea- 
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surement bins (+0"50, +0'/25, O'.'OO, and -0"23) in the Markarian 3 STIS spectra show an 
asymmetry: the observed flux of the blueshifted velocity component is less than the flux 
of the redshifted component (e.g. Figure 7). This asymmetry is observed for resonance, 
semi-forbidden, and forbidden lines. We do not see evidence in Figure 8 that the extinc- 
tion is commensurately greater for the blueshifted components compared to the redshifted 
components. The extinction corrected redshifted/blueshifted flux ratios are also typically 
greater than unity. The redshifted components may be brighter simply because they have 
higher emissivity (i.e., more gas) than the blueshifted components. 

We estimated the hydrogen column density along the line-of-sight to the narrow-line 
region clouds using our derived color excess values and the relation of Shull & van Steen- 
berg (1985): N#/=5.2xl0 21 cm~ 2 x E(B-V). We found that the column density varies from 
0.35 xlO 21 cm -2 to 1.8 xlO 21 cm" 2 corresponding to our range in calculated E(B-V) of 0.07 
to 0.36. These values are reasonable for an SO galaxy viewed at i=33 (Tuffs et al. 2004), 
the inclination of Markarian 3 (Schmitt & Kinney 2000). 

4. Emission-Line Diagnostics 

4.1. Ionization Diagnostics 

We used emission-line ratios to constrain the ionization mechanism in the NLR. Veilleux 
& Osterbrock (1987) identified line ratios to distinguish between starburst and AGN con- 
tinua. They used lines that are close in wavelength to diminish the effects of reddening. We 
plot the logarithm of the ratio [N II] A6584/Ha A6563 against that of [O III] A5007/H/3 A4861 
in Figure 9. The solid line in each figure separates the starburst parameter space from that 
for AGN (Veilleux & Osterbrock 1987). The dashed line represents the power-law pho- 
toionization model with solar abundance of Ferland & Netzer (1983). The dimensionless 
ionization parameter varies from 10~ 4 in the lower right to 10 -1 ' 5 in the upper left. The 
ionization parameter is defined as U = Q (H)/4-7rr 2 n e c, where Q(H) is the number of ionizing 
photons emitted by the central source per second, r is the distance between the source and 
the cloud, n e is the electron density at the inner face of the cloud, and c is the speed of light. 
If the density is low enough that collisional de-excitation is negligible (<10 3 cm" 3 for most 
of the observed emission-lines), emission-line ratios will depend primarily on U (Ferland & 
Netzer 1983). For example, the [O III] A5007/H/3 A4861 and the [O III] A5007/[O II] A3727 
ratios increase as U increases. 

The observations occupy a small region of parameter space that is approximately 0.2 dex 
higher than the solar-abundance power-law photoionization model of Ferland & Netzer 
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(1983). That the observations occupy a small locus of parameter space indicates that the 
same ionization mechanism, AGN photoionization, applies for all the gas measured along 
the slit. This is consistent with the results of Kraemer & Harrington (1986) and of Sako et 
al. (2000) for Markarian 3. The difference between the model and the data in the logarith- 
mic ratio corresponds to a factor of (0.2 dex=10°' 2 =) 1.6. This may be an indication that 
the N/O abundance ratio is somewhat higher than solar (for solar abundances see Lambert 
(1978)). Such an offset is typical of Sy 2 NLR (see Figure 12.1 in Osterbrock (1989)). 
An enhanced N/O ratio may be an indication of a greater than solar nitrogen abundance 
combined with a less than solar oxygen abundance. Oxygen may be depleted onto grains in 
a dusty nebula. We will investigate this possibility in Paper II. 

We also examined the line ratios log([0 I] A6300/H« A6563) versus log([0 III] A5007/H/3 A4861). 
This indicator occupies a small region of its AGN photoionization parameter space similar 
to that of the [N II] A6584/Ha A6563 diagnostic. It does not show the same offset from the 
Ferland & Netzer (1983) photoionization model however, suggesting that the abundance for 
oxygen is close to solar. This is also true of the plots in Osterbrock (1989). 



4.2. The Ionizing Continuum 

The He II A4686/H/5 A4861 ratio is sensitive to the shape of the ionizing continuum 
(Kraemer 1985; Kraemer et al. 1994) if the NLR clouds are ionization-bounded. The line 
ratio is lower with a steep spectrum than with a flat spectrum because there are less He 4 " 
ionizing photons in the former case. We assumed that neutral hydrogen is the only absorber 
of photons with energies between 13.6 eV and 54.4 eV, and that He + is the only absorber of 
photons with energies greater than 54.4 eV. For a power-law spectrum with spectral index a 
the ratio of the number of hydrogen ionizing photons, Q(H°), to the number of He + ionizing 
photons, Q(He + ), is 

R _ Q(H°) N(He) I (HP) j (4686) ^ fl\ a 

Q(He+) N(H) 1(4686) j(H(3) ~ \A J 1 ' 

where N(He)/N(H) has the assumed solar value of ~0.1, I(H/9)/I(4686) is the extinction- 
corrected measured line ratio, and j (4686) /j (H/3) is the ratio of volume emissivities obtained 
from Osterbrock (1989). We summed the line fluxes measured in all bins for He II A4686 and 
H/3 and found a ratio I(H/9)/I(4686) of 4.57. This yields a=1.36. This agrees with the value 
of a=lA we found for NGC 1068 (Sy 2) (Kraemer & Crenshaw 2000b), NGC 5548 (Sy 1.5) 
(Crenshaw & Kraemer 1999), NGC 3783 (Sy 1) (Kraemer, Crenshaw, & Gabel 2001), and 
NGC 4151 (Sy 1.5) (Kraemer & Crenshaw 2000a) in the spectral range 13.6 eV < E < 1 keV. 
However, if there is a preponderance of matter-bounded clouds, the intrinsic continuum could 
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be steeper. 

4.3. Extinction-Corrected Line Ratios vs. Angular Separation from the 

Kinematic Center 

The ionization-parameter sensitive line ratios can be used to probe the electron density 
as a function of position within the NLR. If n e is constant throughout the NLR, the ioniza- 
tion parameter and its diagnostic line ratios decrease as r~ 2 in the NLR (see the definition 
of the ionization parameter, U, in §4.1). If n e varies with radius, the radial dependence of U 
(and its diagnostic line ratios) will be modified since Uoc(l/r) 2 x(l/n e ). Figure 10 shows the 
extinction corrected flux ratio [O III] A5007/H/3 A4861 vs. projected position in the NLR 
along the STIS longslit aperture. A similar plot for the ratio [O III] A5007/[O II] A3727 is 
shown in Figure 11. Using the bi-cone geometry of Ruiz et al. (2001) (inner and outer 
opening angles of 15° and 25°, respectively, and tilted out of the plane of the sky by 5° 
towards the observer in the east) and of Schmitt & Kinney (1996) (PA. =70°) we can de- 
termine the radial distances of the NLR clouds in Markarian 3. We use the solar abundance 
photoionization model of Ferland & Netzer (1983) to convert the line ratios to photoion- 
ization parameter values. The electron density is a constant 10 3 cm -3 for this model. The 
ionization parameter values derived from the two different line ratios are shown as functions 
of radial distance in Figure 12 and Figure 13. We show curves on the figures representing U 
corresponding to different radial distributions of n e . For the curve labeled "2.0", n e oc(l/r) 2 
(and U is constant with radius), while n e is constant with radius (and Uoc(l/r) 2 for the curve 
labeled "0.0". Note that in the former case, the curve is flat because the radial dependence 
of the electron density cancels the explicit radial dependence of the ionization parameter. 
The values of U at large radius in both Figures 12 and 13 lie between the two plotted curves 
of n e oc(l/r) 2 and of n e =constant. This indicates that the electron density in the NLR 
decreases with radius from the nucleus. The uncertainty in the ionization parameter is large, 
but the overall trends in both figures suggest that the decrease is slower than (1/r) 2 . A 
patchy gas distribution along the slit in the NLR might be responsible for deviations from 
the curves. Nelson et al. (2000) and Kaiser et al. (2000) observed similar, although 
smoother, trends in STIS observations of NGC 4151 and concluded that the electron density 
decreases slower than r~ 2 in its NLR. 

However, collisional de-excitation of the + ions could contribute to the higher line ratio 
close to the nucleus. The critical densities of the upper-levels responsible for the [O II] A3727 
doublet are 1.6xl0 4 cm" 3 for 2 D 3/2 and 3.1xl0 3 cm" 3 for 2 D 3/2 (Osterbrock 1989). These 
are, respectively, factors of 44 and 226 less than the critical density (7.0 xlO 5 cm -3 ) for the 
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upper level ( 1 D 2 ) responsible for the [O III] A5007 line (Osterbrock 1989). If the electron 
density near the center of the NLR were greater than the critical densities for the 2 D 3 / 2 
and 2 D 3 / 2 levels of + and less than the critical density of the 1 D 2 of +2 , then the flux of 
the [O II] A3727 doublet will drop, and the [O III] A5007/[O II] A3727 ratio will increase. 
We find from the density and temperature sensitive emission-line ratios described below, 
however, that the electron densities in the NLR are probably lower than those conducive 
to collisional de-excitation. In either case, with or without collisional de-excitation, the 
drop in the ionization parameter sensitive diagnostics with increasing distance from the 
nucleus indicates a radial dependence on electron density. Photoionization modeling will 
be employed in Paper II to further investigate the physical conditions that produce these 
observed gradients. 

The line ratio [O III] AA (5007+4959) /A4363 is a sensitive temperature diagnostic for a 
photoionized gas in the density range 10 cm -3 < n e < 10 4 cm" 3 . The theoretical value of 
this ratio decreases with increasing temperature. We detect no obvious trend with spatial 
position for this diagnostic. The scatter in the [O III] ratio corresponds to an uncertainty in 
the temperature ranging from 12000 K to 17000 K. A similar temperature range is found in 
the NLR of NGC 4151 (Nelson et al. 2000). 

Although the energy required to produce most of the ions observed in the Markarian 3 
NLR is less than 65 eV, two ions, Fe +6 and Ne +4 require energies of 100 eV and 97 eV, re- 
spectively. The spatial profiles of emission lines from these species exhibit striking differences 
from the lower ionization lines at neighboring wavelengths. The flux drops rapidly from the 
nucleus to the eastern- most detectable component. This is most dramatically illustrated by 
comparing the [Fe VII] A5721 and A6087 lines with their neighboring lines, [N II] A5755 and 
He I A5876, in the mode G750L spectrum (Figures 5 and 6). The fact that lines close in 
wavelength do not show the same spatial morphology suggests that this is not a reddening 
effect. The flux measurements for the [Fe VII] A5721 line relative to ¥L(3 A4861 are plotted 
in Figure 14. Figure 15 shows the flux measurements of [Ne V] A3426 relative to H/3 A4861. 
The asymmetric spatial profile of the emission-lines from these high-ionization potential ionic 
species suggests that an absorber located between the BLR and the NLR blocks the contin- 
uum above ~65 eV that would be required to ionize Fe +5 and Ne +3 in the eastern side of 
the NLR. Similar effects have been suggested for NGC 4151 (Alexander et al. 1999) and 
NGC 1068 (Alexander et al. 2000). 
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5. The Observed Continuum 

5.1. Spectral Decomposition of the Continuum 

We examined the observed continuum in the same 1'.'8 region spanned by the bins defined 
for the emission-line measurements. Since we found a reddening gradient along the STIS slit 
increasing from west to east from emission-line diagnostics, we surmised that the intrinsic 
continuum may also be similarly reddened. We made separate continuum extractions for the 
east and for the west. The western extraction is the sum of the detector rows corresponding 
to the bins labeled -0'.'8, -0'.'5, -0'.'2, and O'.'O in Tables 2 and 3. The eastern extraction 
is the sum of the rows corresponding to bins 0'.'3, 0'.'5, and 0'.'8. For both extractions we 
exclude those regions containing emission lines in the row sum. After applying a foreground 
reddening corresponding to E(B-V)=0.08 (LMC extinction curve [Koornneef & Code 1981]) 
to the western extraction, we found that its shape matched that of the eastern extraction. 
The flux level of the eastern extraction is 1.4 times greater than that of the extinction- 
transformed western extraction over the entire observed wavelength range. The eastern 
spectral extraction and the extinction-transformed western extraction are shown with the 
model components for the eastern extraction in Figures 16-17. Note that the extinction- 
transformed western extraction is scaled by 1.4 for comparison with the eastern extraction. 

We modeled the observed continuum as a combination of host galaxy stellar light, scat- 
tered light from the active nucleus, and recombination emission from H + and He +2 in the 
NLR. Stellar light solely from an old population has insufficient far-UV (STIS mode G140L) 
flux to match the observed spectrum. We determined the relative contribution of the stellar- 
light and scattered AGN-light components to our observed spectrum using the Kinney et al. 
(1996) spectral templates for normal galaxies and a simple power-law spectrum. We used 
the Kinney et al. (1996) SO galaxy template for the normal galaxy component. We assumed 
that the spectral index of the power-law component in the observed spectral region is a—1. 
The recombination continuum was determined using the recombination-line and continuous- 
emission calculations described in Osterbrock (1989). The input variables for these calcu- 
lations were the sums of the extinction corrected Hf3 A4861 and He II A4686 measurements 
from Table 3 in those bins corresponding to eastern (+0'. / 3 to 0'.'8) and western (-0'.'8 to O'.'O) 
extractions. We applied the Galactic extinction (corresponding to E(B-V)=0.19 [Schlegel et 
al. 1998]) in the direction of Markarian 3 to all components of the model spectrum using 
the Savage & Mathis (1979) extinction curve. We also applied an extinction corresponding 
to E(B-V)=0.08 using the LMC curve (Koornneef & Code 1981) to all components of the 
model spectrum to match the eastern extraction. We normalized the stellar and power-law 
components of the model spectrum so that they matched the data at the 8125 A mea- 
surement bin of the eastern extraction. We scaled and added the two components together, 
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varying the scale factor to determine the ratio of host galaxy to scattered AGN light that 
best matched the overall shape of our observed continuum. The signal-to-noise is too low to 
permit matching absorption lines between the data and the galaxy template. The recombi- 
nation continuum model components were added to the power-law + host galaxy model. A 
host galaxy to power-law light ratio of 3:1 at 8125 A yielded a good match to the data, but 
we found that additional extinction intrinsic to the Markarian 3 host galaxy was required 
to give the best fit at all wavelengths. Applying additional extinction corresponding to E(B- 
V)=0.20 with the LMC curve (Koornneef & Code 1981) to the host galaxy component and 
to the recombination continuum components gave the best fit to the eastern extraction. The 
reddened model components for the best fit to the eastern extraction are shown with the 
eastern and western extractions in Figure 16. The sum of the three model components are 
shown with the data in Figure 17. 

The uncertainty in the host galaxy fraction {i ga iaxy = L ga i axy /(Lg alaxy +L scattered _ AGN , 
where L ga i axy and L scatte red- agn are the host galaxy and power-law continua, respectively) 
is shown by the gray shaded area in Figure 17. The galaxy fraction at 8125 A ranges from 
0.65 to 0.85 across this region. We estimated the uncertainty in the host galaxy extinction 
for a given host galaxy template and host galaxy extinction value for a fixed galaxy fraction 
at 8125 A of 0.75. Varying the host galaxy color-excess from ±0.05 matches the scatter 
in the smallest error bars at wavelengths longer than 4000 A. A variation of ±0.10 matches 
the scatter in the largest error bars. Variations at these levels have no effect on the model 
spectrum at wavelengths less than 4000 A, where the scattered continuum from the hidden 
AGN dominates. 

Our results for continuum decomposition have a lower galaxy fraction than the previ- 
ous studies of Koski (1978), Malkan & Oke (1983), and Tran (1995a). The wavelength 
dependence of the galaxy fraction on our continuum model normalized to 0.75 at 8125 A is 
shown in Figure 18. We estimated that the host galaxy fraction reaches 0.8 at the longest 
wavelengths (>9500 A) in our model and drops to 0.65 near 4500 A. The fraction drops 
off sharply towards shorter wavelengths, amounting to only 0.05 at 2600 A. The fraction 
gently decreases to 0.02 at far-UV wavelengths. Koski (1978) reported a galaxy fraction at 
4861 A of 0.78. Malkan & Oke (1983) found a galaxy fraction of 0.88 at 5175 A and Tran 
(1995a) computed a fraction of 0.88 at 5500 A. We found galaxy fractions close to 0.60 at 
4861 A and 5171 A, and 0.70 at 5500 A. These differences are likely due to differences in 
aperture size. The previous studies used apertures of width 2'/4, 1", and 2'.'4, respectively, 
which admit more host galaxy light than the STIS C/.'l wide slit. If the same amount of 
scattered nuclear continuum light is observed through all apertures, the host galaxy fraction 
derived from a large aperture observation will be higher than that derived from a small 
aperture observation. 
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We found that the total flux of our observed continuum, integrated in wavelength from 
1164 A to 9670 A, is 2x10 12 ergs s 1 cm 2 . When corrected for Galactic extinction only, 
this value is 3.5xl0™ 12 ergs s" 1 cm -2 . At a distance of 53 Mpc, this total luminosity is 
l.lxlO 42 ergs s" 1 . 

The scattered AGN continuum luminosity can be derived from our power-law fit by 
recalling that the western component has an unreddened power-law shape with a—1 and 
the eastern component is a reddened (E(B-V)=0.08) version of the western component and 
is 1.4 times brighter. When we add these two components together, and integrate from 
1164 A to 9670 A we find that the luminosity is 5.8xl0 41 ergs s _1 . The total luminosity 
of the unreddened AGN continuum scattered light integrated from 1164 A to 9670 A is 
7.4 xlO 41 ergs s _1 . 



5.2. The Covering Factor 

We estimated the covering factor, C, of the NLR gas using the relationship between H/3 
luminosity and ionizing continuum given by Netzer & Laor (1993) 

L C 

L(Hf3) ~ 7.4 x 10 13 -^ ergs s' 1 (3) 

where L\ is the monochromatic luminosity at 1 Ryd in ergs s _1 Hz" 1 and a is the power-law 
index above 1 Ryd. This equation applies to an optically thick gas under case B recombi- 
nation conditions. We summed our extinction-corrected H/3 A4861 flux measurements listed 
in Table 3 then scaled by the distance factor AttD 2 to find a total H/3 A4861 luminosity of 
7.2±0.4xl0 40 ergs s -1 . To find the luminosity at 1 Ryd, we assumed that the ionizing con- 
tinuum has the form of a two component power-law with a break-point at u=(0.2 keV)/h, 
where h is Planck's constant. We used the spectral index of the ionizing continuum de- 
rived in §4.3 of a— 1.36 for the range (13.6 eV/h)< v <(0.2 keV/h). Turner et al. (1997) 
used the ASCA spectrum of Markarian 3 to estimate that the unabsorbed, intrinsic X-ray 
spectrum has a spectral index a—1 over the range (0.2 keV/h)< v <(50 keV/h) and that 
its integrated flux from 2 keV to 10 keV is ~10 44 ergs s" 1 . Using this normalization, we 
found that Lx=5.0xl0 28 ergs s _1 Hz™ 1 , and the covering factor is 2.7±0.15%. Netzer & Laor 
(1993) found typical NLR covering factors for samples of Seyfert 1 galaxies and radio- loud 
quasars to lie between 1% and 4%. We note that the total ionizing continuum luminosity 
from 13.6 eV to 50 keV using this two-component model is 6.2 xlO 44 ergs s _1 . 

The covering factor value derived from the R{3 A4861 measurement should be a lower 
limit since our equivalent aperture of size 0f!l x l'/8 does not sample the total emitted 
B.p A4861 flux. From purely geometrical considerations, the maximum possible covering 
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factor for our adopted hollow bi-cone geometry (with opening half angles of 15° and 25°) 
is (Af2/47r) xl00% = 6%. At our adopted distance to Markarian 3 the slit width is 26 pc 
and our H/5 measurement region covers ~30% of the NLR bi-cone. The upper limit to the 
covering factor constrained by the superposition of the STIS slit on our bi-cone model is 6% 
x 30% = 2%. This is on the order of our photon-counting derivation lower-limit. Based on 
this analysis the covering factor within the slit is unity given the bi-cone geometry. However, 
if we have underestimated the value of L 1; the covering factor may actually be lower (see 
below) . 

We estimated the fraction of intrinsic non-ionizing AGN continuum emission scattered 
into our line-of-sight. We extrapolated our power-law fit of the observed continuum (cx—1) 
to z/=(13.6 eV)/h and scaled it to the value of L x estimated for the ionizing continuum. The 
inverse of this scale factor, 0.2%, is the scattering fraction. 

We compared the scattering fraction of the continuum with that for BLR emission. 
We detected scattered C IV AA1548,1551 broad-line emission in a large bin (l'/8) extraction 
equivalent to the full spatial range of the seven emission-line measurement bins. These 
broad-line wings are indicated by the heavy gray line in Figure 2. Other permitted lines in 
the STIS spectrum may have scattered broad components, but they also have neighboring 
lines which make such identification difficult. Of the permitted lines observed in the STIS 
spectrum, C IV AA1548,1551 and Mg II AA2796,2803 are the most isolated. We do not 
observe scattered broad-line flux from Mg II. 

We followed the procedure of §3.1 to measure the fluxes of the various components 
of the blended C IV AA1548,1551 lines. We assumed that each narrow line consists of two 
kinematic components: one redshifted and the other blueshifted. We identified the individual 
narrow line widths and velocities, relative to the systemic velocity, by using the [O III] A5007 
line as a template. We fitted the remaining flux above the continuum with a Gaussian 
function. The scattered broad-line component flux is 1.25xl0 -14 ergs s _1 cm -2 . The broad- 
line FWHM of 24 A corresponds to a velocity dispersion of ±2200 km/s. After correcting 
for Galactic extinction (Savage & Mathis 1979) (corresponding to E(B-V)=0.19, [Schlegel 
et al. (1998)]) we found a flux of 5.2xl0~ 14 ergs s _1 cm -2 . Accounting for the 53 Mpc 
distance to Markarian 3 the luminosity of the scattered line emission is 1.6xl0 40 ergs s _1 . 
We note that the equivalent width of the broad-line component is ~125A. This is consistent 
with the measurements of Wu, Boggess & Gull (1983) for Sy 1 galaxies. 

We estimated the intrinsic C IV AA1548,1551 luminosity by assuming a fixed ratio of 
hard X-ray flux to broad-line flux. We scaled the Markarian 3 hard X-ray flux by that ratio 
determined for Fairall 9, a Sy 1 galaxy with relatively unabsorbed AGN X-ray continuum 
emission. We obtained the hard X-ray flux, f2_io fceV=2xlO~ n ergs s" 1 cm -2 , from the TAR- 
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TARUS database. Wu, Boggess & Gull (1983) listed a C IV AA1548,1551 flux of 7.31 xl(T 12 
ergs s _1 cm -2 for Markarian 3. Scaling our adopted hard X-ray luminosity of 10 44 ergs s^ 1 
for Markarian 3 by the ratio of C IV AA1548,1551 to hard X-ray flux for Fairall 9 yielded an 
estimated C IV AA1548,1551 luminosity of 3.7xl0 43 ergs s _1 for Markarian 3. Dividing our 
measurement of the scattered line luminosity by this estimated intrinsic luminosity indicates 
that the line emission scattering fraction is 0.04%. This is a factor of 5 less than our estimate 
for the continuum scattering fraction. 

We can obtain a better agreement between the continuum and line emission scattering 
fractions if the ionizing continuum has a steeper spectral index for the range (13.6 eV/h)< 
v <(0.2 keV/h). A larger value for L 1 will lower the scattering fraction for the non-ionizing 
continuum emission. We hold the hard X-ray normalization fixed at 10 44 ergs s _1 with a—1. 
We used an X-ray to optical energy index of aox=l-5 (Turner et al. 1997) to predict the 
intrinsic luminosity at 2500 A. This energy index is the modal value of the aox distribution 
for the Braccesi BF quasar subsample studied by Wilkes et al. (1994). We used the non- 
ionizing continuum index of a—1 to extrapolate from 2500 A to 912 A(13.6 eV). We then 
found that the spectral index for the range (13.6 eV/h)< v <(0.2 keV/h) is a ~2. This is 
the same index that Kraemer & Harrington (1986) used to model the Markarian 3 NLR 
emission lines. The new value of L x is 2.78xl0 29 ergs s -1 . Now the scattered fraction of 
continuum emission is 0.04%, which is consistent with the scattered fraction of intrinsic 
line emission. We note that the new values of spectral index and ionizing luminosity at 
1 Ryd yield a covering factor for the NLR gas of >0.7±0.04%. This value is well below the 
geometrical upper limit of 2%. The integrated ionizing luminosity from 13.6 eV to 50 keV 
is now 1.2xl0 45 ergs s _1 . 

The discrepancy between the spectral index derived from the He II A4686/H/3 ratio and 
that derived from extrapolation using the X-ray to optical energy index may indicate that 
the Markarian 3 NLR contains some combination of matter-bounded and ionization-bounded 
clouds. The former derivation applies only to ionization-bounded clouds. Additionally, the 
ionizing continuum might not be well characterized by a simple power law. It may have 
a "big blue bump" in the EUV (Sun & Malkan (1989); Krolik et al. (1991)). We will 
investigate these possibilities using detailed photoionization models in Paper II. 



6. Orientation of the NLR Bi-Cone With Respect to the Host Galaxy 

Schmitt & Kinney (2000) observed an ellipticity (e=l-b/a) of 0.159 in a 60 fim image 
of Markarian 3. This corresponds to an inclination of 33° (i = arccos(b/a)). However, they 
were unable to determine the orientation of the Markarian 3 galactic disk. Galactic disk 
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orientation may be determined using rotation curves derived from neutral hydrogen (21 cm) 
observations, but no such information for Markarian 3 is currently available in the literature. 

Schmitt & Kinney (2000) suggest that dust lanes should lie only on the closer side 
of a galactic disk, as they would be delineated by background bulge light. We searched 
for dust lanes in an archival HST/WFCPC2/F606W image obtained by Malkan, Gorjian, 
& Tarn (1998) (exposure time is 500 seconds). After cleaning the cosmic-rays on the single 
archival image, we created an unsharp- masked image of Markarian 3 to enhance high contrast 
features that are overwhelmed by the bulge light of Markarian 3. The unsharp-masking 
process requires dividing an image by a median filtered copy of the same image. We used 
a 31x31 pixels squared (l'/43x l'/43) two-dimensional median filter. The unsharp-masked 
F606W image is shown in Figure 19. The display in Figure 19 is such that emitting-material 
appears white and obscuring material (i.e., dust lanes) appear black. The horizontal bar 
corresponds to 1 kpc or 3" 8. Note the preponderance of concentric dust lanes found between 
4" and 12" north-east of the nucleus, and the lack of such lanes in the south-west. This is 
consistent with the reddening gradient observed along the STIS slit that increases from west 
to east as shown in Figure 8. Based on the orientation of the galaxy disk major axis (P.A.=28° 
east of north [Schmitt and Kinney 2000]) and using the bulge back-lighting assumption, we 
would expect such arcs to be oriented perpendicular to the minor axis on the near side of 
the galaxy, not perpendicular to the major axis as they appear in Figure 19. However, the 
presence of dust lanes in the east and not in the west and the bulge backlighting hypothesis 
suggest that the host galaxy disk is tilted towards the observer in the east along the STIS 
slit (P.A.=70° [Schmitt and Kinney 1996]). 

We used these clues to the host galaxy inclination and the results of our line diagnostic 
and continuum studies to develop a physical picture of the Markarian 3 NLR. The reddening 
gradient derived from the He II AA1640,4686 line emission shown in Figure 8 suggests the 
presence of a foreground screen between the observer and the NLR, with increasing extinction 
from west to east. An extinction gradient is supported by the continuum observations as 
well: the continuum shape of the western extraction matches the eastern extraction only 
after applying an amount of extinction corresponding to E(B-V)=0.08 using the Koornneef 
& Code (1981) LMC curve. We infer that the foreground screen is dust within the inclined 
disk of the host galaxy and argue that the plane of the host galaxy is tilted towards us in 
the east. We adopt the disk inclination of 33° from Schmitt & Kinney (2000). Ruiz et al. 
(2001) showed that the NLR bi-cone is tilted towards us in the east as well, although at a 
smaller angle (5°). The eastern bi-cone is more heavily reddened because we view it through 
the dust within the galactic plane. The western bi-cone is less reddened since it lies above 
the plane of the galaxy, or at least the line-of-sight intersects less planar dust than that 
towards the eastern bi-cone. This configuration is consistent with the findings of Elmegreen 
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& Block (1999). They suggest that the near side of an inclined disk galaxy will be redder 
than the far side because the light source is behind most of the disk dust. 

In addition to reddening in the east, the contrast in the extinction gradient might be 
enhanced by preferentially scattered blue continuum emission in the west. If the host galaxy 
disk is tilted away from us in the west, dust within the disk and throughout or behind the 
NLR might scatter blue continuum emission into our line-of-sight. Kishimoto et al. (2002) 
found evidence of more scattering (possibly due to dust) in the west than in in the east from 
a polarized-flux color-map (HST/WFPC2 F275W-F342W) of the Markarian 3 NLR. 

Figures 20 and 21 are schematic diagrams of our physical picture of the Markarian 3 
NLR shown from two different viewing angles. We show the relative orientation of the 
components, not their relative spatial scales. The host galaxy disk position-angle measured 
east from north is 28° and its inclination is 33° (Schmitt & Kinney 2000). The host galaxy 
disk is tilted towards the observer in the east. The bi-cone orientation parameters from Ruiz 
et al. (2001) are: position-angle (measured east from north) = 72°, outer-edge opening 
half-angle = 25°, and inclination = 5° out of the plane of the sky towards the observer in 
the east. These parameters are derived from the [O III] A5007 morphology and kinematics 
of the inner 2-3 " of the NLR. In Figure 20 the point-of-view is that of the observer such 
that the plane of the page is parallel to the plane of the sky. North is at the top of the figure 
and east is to the left. In Figure 21. The vertical line is parallel to the observer's projected 
line-of-sight. The observer views the system from the top of the figure. We assume from 
Figure 8 of Ruiz et al. (2001) that most of the NLR gas is confined within the inner (15° 
[Ruiz et al. 2001] not shown in the figures) and outer half-opening angles of the bi-cones. 

Ruiz et al. (2001) found that there were more redshifted points in the East further out 
in the ENLR than in the West at comparable distances from the nucleus (see their Figure 4). 
If the motions of the ENLR are governed by the galactic disk rotation, these observations and 
the inference that the disk is inclined towards the observer in the East would be consistent 
with a clockwise (north-to-west) disk rotation as viewed by the observer. Kotilainen & Ward 
(1997) presented color maps of Markarian 3. Their V-R color map shows a red spiral feature 
extending approximately ±4'/5 from the nucleus. If this feature has a kinematic origin and if 
the spiral arms are trailing the rotation of the galaxy, it would support the hypothesis that 
the galaxy rotation runs clockwise from north to west. 
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7. Summary 

We used HST/STIS low- resolution longslit spectra to study the physical conditions of 
the Markarian 3 NLR. We found an extinction gradient in the emission-line and continuum 
flux increasing from west to east along the STIS slit. The extinction is best characterized 
with an LMC-type Koornneef & Code (1981) curve. We developed a physical picture of the 
NLR and host galaxy geometry using the extinction gradient, the NLR-gas kinematic model 
of Ruiz et al. (2001), and the host galaxy inclination of Schmitt & Kinney (2000). We infer 
that the host galaxy disk is tilted towards us in the east, providing a dustier line-of-sight to 
the eastern NLR bi-cone than that to the western bi-cone. 

Using emission-line diagnostics for the dimensionless ionization parameter (U) we found 
evidence for a large scale gradient in electron density in the NLR. The density appears to 
decrease with radius, although slower than (1/r) 2 . Small scale deviations from the trend in- 
dicate that the gas distribution within the NLR may be patchy. Diagnostics for temperature 
show no trend with position along the STIS slit. The scatter in the diagnostic suggests an 
uncertainty in temperature ranging from 12,000 K to 17,000 K. Photoionization diagnostics 
are consistent with AGN power-law ionizing radiation and show no evidence for a starburst 
in Markarian 3. We found evidence that the NLR nitrogen to oxygen abundance ratio is a 
factor of 1.6 above the solar value. 

We decomposed the observed continuum as SO host galaxy stellar light, scattered AGN 
power-law continuum (a=l), and recombination continua from H + and He +2 . The ratio of 
host galaxy continuum to scattered power-law continuum is 3:1 at 8125 A in our aperture. 
We derived the ionizing continuum (E>13.6 eV) as a two-step power law with a breakpoint at 
0.2 keV. We used the Turner et al. (1997) estimate of the unabsorbed X-ray spectrum above 
the breakpoint, where a=l and the 2-10 keV integrated luminosity is 10 44 ergs s -1 . We found 
a=2 for the spectral range 13.6 eV < E < 0.2 keV. When we extrapolated from the ionizing 
continuum to lower energies (E < 13.6 eV) and compared the result with the scattered- 
AGN component of our observed continuum fit we found that the observed scattered-AGN 
continuum in the UV is 0.04% of the intrinsic AGN continuum. This is consistent with 
our estimate for the percentage of intrinsic C IV AA1548,1551 line emission scattered into 
our line-of-sight. Given our adopted bi-cone geometry, the STIS aperture through which we 
observed the NLR (52 // x0'. / l), and the ionizing continuum slope a=2 at 13.6 eV, the covering 
factor of the NLR gas is 0.7%<C<2%. 

In Paper II we will use the physical parameters derived from the emission line diagnostics 
and the inferred ionizing continuum as input to the photoionization model code CLOUDY 
(Ferland 1988) to confirm our physical picture of the NLR. These models will put tighter 
constraints on sizes, densities, and spatial distribution of gas clouds in the NLR. 
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N 




Fig. 1. — STIS target acquisition image obtained through the 
F28X50LP filter with 52"x0'. / l slit superposed. The horizontal bar 
drawn below the NLR is 1" or 257 pc in length. The slit is oriented 71° 
east of north. East is towards the bottom of the figure and North is 
towards the left. The arrow indicates the NLR kinematic center derived 
by Ruiz et al. (2001) and the position of the putative hidden nucleus 
identified by Kishimoto et al. (2002). 
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Fig. 2. — l'/8 extraction from G140L spectral image. The heavy gray line at 
C IV AA1548,1551 shows the fit to the broad-line flux scattered into our line-of-sight 
that is discussed in §5.2. The l'/8 spectral image extraction is displayed at the bottom 
of the panel. 
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Fig. 3. — l'/8 extraction from G230L spectral image. 
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Fig. 4. — l'/8 extraction from G430L spectral image. 
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Fig. 5. — l'/8 extraction from G750L spectral image. 
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Fig. 6. — Measurement bin diagram shown for 50 emission lines. The spatial orientation 
is west at the top and east at the bottom. The spatial range from the bottom of the 
eastern-most bin to the top of the western- most bin is l'/8. The spectral range shown 
for each line is 15 A for mode G140L (top row), 40 A for mode G230L (second row), 
68 A for mode G430L (third row), and 123 A for mode G750L (bottom row) The target 
acquisition image (see Figure 1) is duplicated at the ends of each row for reference. At 
the right end of each row, the 52"x0'. / l slit boundaries are drawn on the target acquisition 
image. The asterisks next to some of the line labels indicate that more than one emission 
line is present in the box. The arrow indicates the bin defined as O'.'O in this paper. 
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Fig. 7. — Two component Gaussian fit to the [O III] A5007 line in the central bin (position 
= O'.'O). The dashed vertical line indicates the wavelength of [O III] A5007 at the systemic 
velocity of Markarian 3. The dotted line near the bottom of the emission-line shows the 
fit to the continuum. The dashed lines show the fits to the two emission-line components. 
The diamonds show the sum of the emission-line and continuum fits. 
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Fig. 8. — E(B-V) color excess as a function of angular separation from the nucleus 
along the slit. Vertical bars show the uncertainty in E(B-V), and horizontal bars show 
the angular size of the measured bin. Open squares represent redshifted components 
and filled squares represent blueshifted components. There is a clear extinction gradient 
increasing from west to east in the NLR. The zero-point of the position axis for Figures 8- 
15 is the kinematic center determined by Ruiz et al. (2001). 
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Fig. 9. — Diagnostic ratio plot to test the nature of the ionizing continuum. The solid 
line separates starbursts from AGN. The dashed line is the power-law photoionization 
model for solar abundance of Ferland & Netzer (1983). 
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Fig. 10. — The ionization parameter sensitive line ratio [O III] A5007/H/5 A4861 vs. 
angular separation from the nucleus along the slit. Open squares represent redshifted 
components and filled squares represent blueshifted components. Generally, the ratio 
decreases with distance from the nucleus for both redshifted and blueshifted components. 
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Fig. 11. — The ionization parameter sensitive line ratio [O III] A5007/[O II] A3727 vs. 
angular separation from the nucleus along the slit. Open squares represent redshifted 
components and filled squares represent blueshifted components. The ratio decreases 
with distance from the nucleus for both redshifted and blueshifted components. 
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Fig. 12. — The ionization parameter sensitive line ratio [O III] A5007/H/5 A4861 vs. radial 
distance from the nucleus towards the east and west along the slit. Open squares represent 
redshifted components and filled squares represent blueshifted components. The curves 
represent theoretical ratios for different radial dependencies on electron density. The label 
above each curve refers to the exponential dependence on radius (for example, the curve 
labeled 2.0 represents a line ratio for which n e oc (1/r) 2 ). All curves are normalized to 
a ratio of 0.0013 at a radial distance of 100 pc west. 
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Fig. 13. — The ionization parameter sensitive line ratio [0 III] A5007/[O II] A3727 vs. 
angular radial distance from the nucleus towards the east and west along the slit. Symbols 
and curves are the same as those described in Figure 12. All curves are normalized to a 
ratio of 0.004 at a radial distance of 100 pc west. 
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Fig. 14. — [Fe VII] A5721/H/3 A4861 plotted as a function of angular separation from the 
nucleus along the slit. Note the steep drop in the ratio to the east. 
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Fig. 15. — [Ne V] A3426/H/3 A4861 plotted as a function of angular separation from the 
nucleus along the slit. Note the steep drop in the ratio to the east. 
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Fig. 16. — The eastern and western extracted continuous spectra (crosses) and individual 
model components. Crosses denote average continuum values measured between emission 
lines: the horizontal bar is the width of the measurement region in wavelength, and 
the vertical bar is the ±la uncertainty in the computed average. Eastern and western 
extractions are overplotted to emphasize their similarity after accounting for extinction 
and intrinsic luminosity. The gray crosses represent the eastern extraction corresponding 
to the region encompassed by bins +0'/80 to +0'/30 in Figures 8-15. The black crosses 
represent the western extraction (H-C/OO to -l'/OO) transformed to match the eastern 
extraction using (1) reddening that corresponds to to E(B-V)=0.08 with the Koornneef 
& Code (1981) LMC extinction curve and (2) a multiplicative scale factor of 1.4. The 
dashed line is the power-law continuum («=1) reddened by an amount corresponding 
to E(B-V)=0.08 with the Koornneef & Code (1981) LMC extinction curve and E(B- 
V)=0.19 with the Galactic extinction curve of Savage & Mathis (1979). The dotted line 
represents the SO galaxy template of Kinney et al. (1996) reddened by E(B-V)=0.28 
with the Koornneef & Code (1981) LMC extinction curve and by E(B-V)=0.19 with the 
Galactic curve (Savage & Mathis 1979). The solid gray line is the NLR gas recombination 
spectrum reddened in the same manner as the SO template spectrum. The power-law 
and SO galaxy template spectra were scaled so that f ga iaxy = Q-75 at 8125A and also so 
that the sum of these components with the recombination spectrum matches the data. 
A logarithmic wavelength scale, adjusted to the rest-frame of Markarian 3 is used to 
separate the ultraviolet measurements. 
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Fig. 17. — The eastern and western extracted continuous spectra (crosses) and the best fit 
model (thick gray line). The extracted spectra are identical to those shown in Figure 16. 
The thick gray line shows the variation in the reddened model spectrum (SO galaxy + 
power-law + H + and He +2 recombination continua, see text and Figure 16 caption) when 
the galaxy fraction (i galaxy = L galaxy /(L galaxy +L scattered _ AGN ) is varied between 0.65 and 
0.85 at 8125 A. 
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Fig. 18. — The galaxy fraction (i ga i axy ), or the ratio of host galaxy continuum (Kinney 
et al. (1996) SO template) to the sum of the host galaxy continuum and power- law 
continuum (cx—1) as a function of wavelength for the best fit to the observed continuum. 

^galaxy = ^galaxy / (^galaxy + ^ scattered- AG n) = 0.75 at 8125 A. 
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Fig. 19. — Unsharp- masked version of archival HST/WFPC2/F606W image. The 
unsharp-mask filter size is l'/43. The horizontal bar corresponds to 3'/8 or 1 kpc. The 
display scaling is such that emitting-material appears white and obscuring materials (i.e., 
dust lanes) appear black. Note the presence of dust lanes towards the north-east of the 
NLR and lack of dust lanes to the south-west. The emission in the central region is 
primarily due to emission-lines from the NLR within the F606W bandpass. 
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Fig. 20. — Schematic of the Markarian 3 NLR and its relative orientation to the host 
galaxy disk from the observer's point-of-view. North is up and East is to the left. 
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Fig. 21. — Schematic of the Markarian 3 NLR and its relative orientation to the host 
galaxy disk, which is edge-on in this figure. The projected line-of-sight is from the top 
of the figure and parallel to the vertical bar. East is to the left of the figure. 



Table 1. STIS Long Slit Observations of Markarian 3 a 
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Detector 
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Dispersion 
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F28X50LP 


3 c 
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52X0.1 


1080 


O5KS01020 


CCD 


G750L 
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0.60 
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1.58 


2.35 


0.025 


52X0.1 


3238 



a Observations were obtained under HST proposal identification number 8480 by Kraemer and Crenshaw. 

b Spectral resolution is full-width at half-maximum (FWHM) in A and is related to a (for a Gaussian distribution) by a 
= FWHM / 2.354. 

c Exposure time listed for O5K5S01J6Q is that for the third exposure in the data set. The same exposure is shown in 
Figure 1. 



Table 2. Observed Emission-line Fluxes Relative to B.(3 
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( 0.14) 


( 0.10) 



a Projected angular distance from nucleus in arcseconds. 

b Measurement is blue-shifted (b) or red-shifted(r) with respect to the Markarian 3 rest-frame. 
c Angular width of measurement bin in arcseconds. 

d The la measurement uncertainty is listed in parenthesis below each measurement 
c < indicates that the recorded value is the la upper limit. 
f H/3 fluxes arc in units of 10~ 15 ergs sec -1 cm~ 2 . 
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( 0.02) 


( 0.10) 


( 0.05) 


( 0.07) 


( 0.10) 


( 0.12) 


( 0-14) 


( 0.13) 


( 0.14) 


[SIII] 


9532 


1.28 


1.22 


1.24 


1.12 


1.85 


1.48 


1.48 


1.09 


2.13 


1.03 


1.21 






( 0-12) 


( 0.04) 


( 0.04) 


( 0.19) 


( 0.13) 


( 0.13) 


( 0.22) 


( 0.28) 


( 0.37) 


( 0.24) 


( 0.34) 


w 


4861 


4.55 


30.28 


20.81 


9.03 


14.10 


11.71 


33.79 


18.87 


28.97 


16.77 


38.12 



Table 3 — Continued 



West East 



(1) (2) 



-0.8 a 


-0.5 


-0.2 


-0.2 


0.0 


0.0 


0.3 


0.3 


0.5 


0.5 


0.8 


r b 


r 


r 


b 


r 


b 


r 


b 


r 


b 


r 


0.3 C 


0.3 


0.2 


0.2 


0.3 


0.3 


0.3 


0.3 


0.3 


0.3 


0.3 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 



( 0.38) ( 0.69) ( 0.57) ( 1.41) ( 0.94) ( 0.96) ( 4.56) ( 3.67) ( 4.71) ( 3.59) ( 9.85) 



a Projected angular distance from nucleus in arcseconds. 

b Measurement is blue-shifted (b) or red-shifted(r) with respect to the Markarian 3 rest-frame. 
c Angular width of measurement bin in arcseconds. 

d The la measurement uncertainty is listed in parenthesis below each measurement 
c < indicates that the recorded value is the la upper limit. 
f H/3 fluxes are in units of 10~ 15 ergs sec -1 cm~ 2 . 



